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Abstract 

This paper is devoted to the complete classification of space curves under affine 
transformations in the view of Cartan's theorem. Spivak has introduced the method 
but has not found the invariants. Furthermore, for the first time, we propound a 
necessary and sufficient condition for the invariants. Then, we study the shapes of 
space curves with constant curvatures in detail and suggest their applications in 
physics, computer vision and image processing. 
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1 Introduction 

Classification of curves has a significant place in geometry, physics, mechanics, com- 
puter vision and image processing. In geometrical sense, a plane curve with constant 
curvature, up to special affine transformations may be either an ellipse, a parabola 
or a hyperbola [14J. This classification will be obtained by the concept of invariants. 
Geometry of curves in spaces with dimension > 3 has studied with geometers such as 
Guggenheimer [5], Spivak [H] and etc. The aim was finding the invariants of curves 
under transformations. On the other hand, in 14J, study of space curves in the view 
of Cartan's theorem was started but has not completed yet. 

This paper can be viewed as a continuation of the work ^T^, where the authors began 
the classification of space curves up to special affine transformations. We determine all 
of differential invariants and our method is different from the method of Guggenheimer 
^Corresponding Author. E-mail Address: m_nadjafikhah@iust.ac.ir 
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and other existing methods. Also, for the first time, we prove a necessary and sufficient 
condition for the invariants in order that complete the classification. Moreover, we 
classify the shapes of space curves of constant curvatures which has a wide variety of 
applications in physics, computer vision and image processing. The general form of 
these shapes are exist in [5j, but here we try to discuss them in more details. 

In physics, classification of curves up to affine transformations has a special position 
in the study of rigid motions. Suppose we have a particle moving in 3D space and that 
we want to describe the trajectory of this particle. Especially, each curve in a three 
dimensional space could be imagined as a trajectory of a particle with a specified mass in 
the view of an observer. By classification of curves we can, in fact, obtain conservation 
laws. 

Computer vision deals with image understanding at various levels. At the low level, 
it addresses issues such us planar shape recognition and analysis. Some results on dif- 
ferential invariants associated to space curves are relevant to space object recognition 
under different views and partial occlusion. The evolution of space shapes under curva- 
ture controlled diffusion have applications in geometric shape decomposition, smooth- 
ing, and analysis, as well as in other image processing applications (see, e.g. [8l [9]) 
and similar to recent results for planer shapes. For instance, there are some important 
applications of moving frames method in use of the differential invariant signatures |12j . 
In [Tl[2] there exist some applications to the problem of object recognition and symme- 
try. Also, joint differential invariants has been proposed as noise-resistant alternatives 
to differential invariant signatures in computer vision [3j. Practical applications of the 
derived shapes in the latest section are related to invariant signatures, object recogni- 
tions, and symmetry of 3D shapes via the generalization of them from 2D shapes to 
3D ones. 

In the next section, we state some preliminaries about Maurer-Cartan forms and a 
way of classification of maps with the notable role of Maurer-Cartan forms and Cartan's 
theorem. In section three, classification of space curves in under the action of 
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affine transformations is discussed. Finally, in the last section, we study the shapes of 
space curves with constant curvatures and propose some applications of these shapes 
in physics, computer vision and image processing. 

2 Maurer-Cartan form 

Let G C GL(n, R) be a matrix Lie group with Lie algebra C and P : G ^ Mat(n x n) 
be a matrix- valued function which embeds G into Mat(n x n), the vector space of 
n X n matrices with real entries. Its differential is dPs ■ TbG — > Tp(^)Mat(n x n) ~ 
Mat(n X n). 

Definition 2.1 The 1-form ujb = {P{B)}^^ • cIPb of G is called the Maurer-Cartan 
form. It is often written lo = P^^ ■ dP. The Maurer-Cartan form is in fact the unique 
left invariant £— valued 1-form on G such that tuid ■ Ti^G C is the identity map. 
The Maurer-Cartan form to satisfies in Maurer-Cartan equation duj = —uo A uo. The 
Maurer-Cartan form is the key to classifying maps into homogeneous spaces of G. This 
process needs to the following theorem (for a proof we refer the reader to [B]): 

Theorem 2.2 (Cartan) Let G he a matrix Lie group with Lie algebra C and Maurer- 
Cartan form uj. Let M be a manifold on which there exists a C— valued 1-form (p 
satisfying d<j) = —(f) A (p. Then for any point x £ M there exist a neighborhood U of x 
and a map f : U ^ G such that f* u = (f). Moreover, Any two such maps /i, /2 must 
satisfy fi = Lb ° f2 for some fixed B £ G (Lb is the left action of B on G). 

Corollary 2.3 Given maps fi, f2 ■ M ^ G, then fluj = f2 uj, that is, this pull-back 
is invariant, if and only if fi = Lb o f2 for some fixed B £ G. 

In Theorem 2.2, if M is connected and simply-connected, then the desired map 
/ may be extended to all of M [15]. We suppose that G be the special linear group 
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SL(3, R) as a Lie group and we denote its Lie algebra with sl(3, R). This Lie group is 
not simply-connected, so our achievements are local. 

Definition 2.4 An affine transformation of the Euclidean space R^ is the composi- 
tion of a translation in R'^ among with an clement of the general linear group GL(3, R). 
An affine transformation is called special or unimodular, if its matrix part is an element 
of SL(3, R). The group of special affine transformations is the connected coefficient, 
closed subgroup of the Lie group of affine transformations. 

The following section is devoted to the study of the properties of space curves' 
invariants under the action of volume-preserving affine transformations, i.e., the special 
affine group. The number of essential parameters (dimension of the Lie algebra) is 11. 
The natural assumption of differentiability is C^. 

In the next section, by defining the new curve ac instead of a considered regular 
smooth curve c, we will see that the classification of curves in R^ and in the viewpoint 
of Theorem 3.1 is equivalent to the ones in SL(3, R). Thus we find the Maurer-Cartan 
form of SL(3, R) and then its pull-back via the matrix-valued curve ac- In fact, ac s 
play the role of s in Corollary 2.3. This tends to a complete set of invariants of 
ac as 1-forms on R. The derived invariants in a corresponding manner determines 
curves of SL(3, R). Finally in Theorem 3.3, we find that these invariants also provide 
a necessary and sufficient condition for specifying curves in R^ when we supposed the 
action of special affine group SL(3, R) on R^. 

3 Classification of space curves 

In the present section, we achieve the invariants of a space curve up to special affine 
transformations. Prom Theorem 2.2, two curves in R^ are equivalent under special 
affine transformations, if they differ with a left action introduced by an element of 
SL(3, R) and then a translation. 
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Let c : [a, b] — ^ be a curve in three dimensional space which we call the space 
curve, be of class and 

det{c',c",c"')^0, (1) 

for any point of the domain, that is, we assume that c', c" and c"' are linear independent. 
Otherwise, if for example, c'" depends to c' and c" for some interval [a,b], then we 
can simply observe that the curve c will sit in R^, which is not our main topic of 
investigation. Moreover, we can assume that det(c', c", c'") > for being avoid writing 
the absolute value in calculations. 

For the curve c, we consider a new curve, namely ac{t) : [a,b] — > SL(3, R), defined 

by 

{det(c',c",c'")}^/=^ 

which is well defined on the domain of c into the special linear group SL(3,R). We 
can study the new curve in respect to special affine transformations, i.e. the action 
of special affine transformations on first, second and third differentiations of c. If we 
assume that A is a three dimensional special affine transformation, then we have the 
unique representation A = toB which B is an element of SL(3, R) and r is a translation 
in R^. If two curves c and c be the same up to axi A, c = Ao c, then we have 

-/ O ' -II T3 II ->ll TD III 

c = B o c , c = B o c , c = B o c . 

Also from det B = 1 we obtain 

det{B o c', S o c", 5 o c'") = det{B o (c', c", c'")) 
det(c',c",c'"). 

and so we conclude that ac{t) = B o ac{t) and Uc = Lb o Uc, where Lb is the left 
translation for B G SL(3,R). 



det(c',c",c'") = 
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This condition is also necessary because when c and c are two space curves in which 
ac = Lb o oic for an element B G SL(3, R), then we can write 

ae = {det(c',c",c'")}~'/'(c',c",c'") 

= {det(S o (c', c", c'"))}~'/^ B o (c', c", c'") 
= {det(c', c", c'")}-'/' B o (c', c", c'"). 

Thus c' = B o c' and there is a translation r such that A = r o S and c = Ao c where 
^ is a three dimensional affine transformation. Therefore we have 

Theorem 3.1. Let c and c are two space curves, c and c are the same with respect 
to special affine transformations, i.e. c = Aoc when A = to B for translation r in 
and B G SL(3, R) if and only if Uc = Lb o ac where Lb is a left translation generated 
byB. 

Prom Cartan's theorem, a necessary and sufficient condition foi = Lb o ac {B G 
SL(3,R)) is that for any left invariant 1-form a;* on SL(3, R) we have a^{io'^) = Q;*(a;*). 
It is equivalent to a^(a;) = a* (a;) for natural sl(3, R)-valued 1-form lo = .dP for 
matrix- valued function P which embeds SL(3, R) into Mat3 x 3, the vector space of 
3x3 matrices with real entries, and lo is the Maurer-Cartan form. 

We must compute a*{P^^.dP) which is invariant under special affine transforma- 
tions. Its entries are in fact invariant functions of space curves. It is a 3 x 3 matrix 
form which arrays are multiplications of dt (1-forms on [a, 6]). 

Since q;*(P~^ -dP) = -dac, so we calculate the matrix -a'^ and then multiply 
it by dt to have al{P~^ ■ dP). we have 

a-i = det(c',c",c'")'/'- (c',c",c'")"' 



det(c',c",c'")-'/'- 



/ r>'r>" — r>'r>" r" n'" — r" r'" r" r'" — n" n'" ^ 
C3C2 C2C3 C3C5 C1C3 C2C1 C1C2 
C2C3 — c'3C2 C3c'/c'iC3 d-^d^ — (^^l j 
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which c = {ci,C2, gs)^ as a column matrix be the vector representation of curve c. We 
also have 

[det(c', c", c'")]' = det{c",c", c'") + det(c', c'", c'") + det(c', c",c"") 
= det{c',c",c""). 



Thus we see that 



(4 4' 4"\ 

t-2 ''2 ''2 



V^3 



,// Jll J, 

-0 O'j 



■det(c',c",c'") 



/ J Jt Jtl \ 

Ci Ci Ci 



I // in 

<'2 ''2 ''2 



3 / 



C3 / 



After some computations, finally we find that ^ • a'^ is in the following multiple of dt 



( det (c',c",c"") 
3det(c',c",c"') 







det(c",c'",c"") \ 



dct(i 



V 



1 





det(c',c",c"") dct8c',c'",c"") 
3det(c',c",c"') det(c',c",c"') 

-, 2det(c',c",c"") 

3det(c',c",c"') / 



Clearly, the trace of the last matrix is zero and entries of a*(P~^ • dP) and therefore 
entries of the above matrix, are invariants of the group action. 

Therefore according to Theorem 3.1, two space curves c, c : [a, 6] are the same 

under special afHne transformations if we have 



det(c^c^c^"o 

det(c',c",c'") 

det(c',c",c'") 
det(c',c'",c"") 



det(c',c",c""^ 

II Till tJIII \ 



det(c', c", c 

det(c^^^c_ 
det(c',c",c'") 

det(c',c'",c"") 



det(c',c",c'") det(c',c",c'") ■ 

We may use of a proper parametrization a : [a, b] [0, 1] such that the parameterized 

curve 7 = CO cr"^ satisfies in condition det{'y' {s) , ^" {s) , {w)) = and then entries 
over the principal diagonal of a*{P''^ ■ dP) be zero. But this determinant is in fact 
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the differentiation of det(7'(s), 7"(s), 7'"(s)) and for being zero it is sufficient tfiat we 
assume det(7'(s),7"(s),7"'(s)) = 1. On tfie other hand, we have 

c' = (70 a)' = a' ■ (7' o a) 

c" = {c')'^ ■ (7" o o") + a" ■ (7' o a) 



{a'f ■ (7'" o a) + 3 a' a" ■ (7" o a) + a'" ■ (7' o a). 



Thus we conclude that 



det(c', c",c"') = det{a' ■ (7' o a) , {a'f ■ (7" o a) + a" ■ (7' o a) , 

(o"')^ • (7"' o o") + 3 * a'a" ■ (7" o a) + a'" ■ (7' o a)) 
dct{a' ■ (7' o a) , {a'f ■ (7" o a) , {a'f ■ (7'" o a)) 
= (a')^ •det(7'o(7, 7"o(7 7'"o(7) 

= (-')^ 

The last expression specifies a, namely the special affine arc length, is defined as follows 

ft r -,1/6 



a := det{c' {u) , c" {u) , c" {u}) 

J a 



du. 



So (7 is a natural parameter under the action of special affine transformations, that 
is, when c is parameterized by a then for each special affine transformation A, ^4 o c is 
also parameterized by the same a. Furthermore, every curve parameterized by a up to 
special affine transformations is introduced with the following invariants 



A I./ II III iiii\ J ' ""^ 

Xi = det(c,c ,c ), X2 = det(c,c ,c 



(2) 



We call xi and X2 the first and the second special affine curvatures resp. Thus finally 
we have 



a*Jp-\dP) 



^00 XI ^ 
1 -X2 

1 



da. 
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Theorem 3.2 Every space curve of class satisfying in condition under the 
action of special (unimodular) affine transformations is determined by its natural equa- 
tions xi = Xi('^) o,nd X2 = X2{(^) of the first and second special affine curvatures (0j 
as functions ( invariants ) of the special affine arc length. 

Theorem 3.3 Two space curves c, c : [a, h] of class which satisfy in condition 

([iP are special affine equivalent if and only if Xi = Xi '^''^^ X2 — X2- 

Proof: The first side of the theorem is trivial with respect to above descriptions. For 
the other side, we assume that c, c are curves of class satisfying (resp.) in 

det(c', c", c'") > 0, det(c', c", c'") > 0, (3) 

that is, they are not plane curves. Also we suppose that functions xi aiid X2 are the 
same for these curves. 

By changing the parameter to the natural parameter a (mentioned above), we 
obtain new curves 7 and 7 resp., that determinant ([3]) will be equal to 1. We prove 
that 7 and 7 are special affine equivalent and so there is a special affine transformation 
A sech that 7 = ^4 o 7. Then we have c = Ao c and the proof is complete. 

At first, we replace the curve 7 with 5 := t(7) properly, in which case that S 
intersects 7 where r is a translation defined by translating one point of 7 to one point 
of 7. We correspond to £ [a, b] to the intersection of 6 and 7, thus 6{to) = j{to)- One can 
find a unique element B of the general linear group GL(3, R) such that maps the frame 
{S' (to), 6" {to), 5"' (to)} to the frame {f(to), 7"(to), f "(^o)}- So we have B o 6'{to) = 
^'{to), B o 6" (to) = 7"(to)) and B o 5"' {to) = 7"' (to)- B is also an element of the special 
linear group SL(3,R); since we have 

det(7'(to),7"(to),7"'(io)) = det{5' {to),5" {to),5"' {to)) and 
dei{5'{to),5"{to),5'"{to)) = det(Bo (7'(to),7"(to),f"(to))), 

and thus det(i?) = 1. If we prove that r] := Bo5 \s equal to 7 on [a, 6], then by choosing 
A = T o B, there will remain nothing for proof. 
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For curves r] and 7 we have (resp. 



^00 X? ^ 



/ // /// \ / 



(in 



I I II iii\ 
[r] ,7] ,7] 



1 
1 



V 

X2 



and 



(f,7",7'")' = (7',7",f") 



x! ^ 



1 
1 



7 

X2 



J 



Since xi ^i-nd X2 remain unchanged under special affine transformations, so we have 
Xi = Xi = Xi and X2 — xl — X2' therefore, we conclude that t] and 7 are solutions of 
ordinary differential equation Y"" + X2 Y" ~ Xi = where Y depends to t. Because 
of the same initial conditions 

r]{to) = B o 5{to) = 7(io), ??'(to) = B o 5' {Iq) = f(to), 
r/"(to) = B o 5" (to) = 7" (to), V"{t) = B o r (to) = 7"' (to), 



and the generalization of the existence and uniqueness theorem of solutions, we have 
r/ = 7 in a neighborhood of to that can be extended to the whole [a, h]. ^ 



Corollary 3.4 The number of invariants of special affine transformation group acting 
on is two which is the same with another results provided by other methods such as 

m- 

The generalization of the affine classification of curves in an arbitrary finite dimen- 
sional space has been discussed in |10j and a complete set of invariants with a necessary 
and sufficient condition of them for classifying curves up to affine transformations has 
been derived. 
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4 Geometric interpretations applied to physics and com- 
puter vision 



In the present section, the geometric interpretation of the first and the second special 
affine curvatures and their apphcations in physics and computer vision is discussed. 
Since every curve parameterized with special affine arc length a and with constant first 
and second affine curvatures xi X2 fulfilled in relation a'^{cr) = ac{cr).{b), for some 
b £ sl(3, R) via the right action of the Lie algebra. In fact, we assumed that the action 
of the Lie group be the left action [11]. Whereof, Maurer-Cartan matrix of SL(3,R) is 

XI 

1 -X2 



( 



a base for Lie algebra sl(3,R) and one can write (^cW) — ctc(c)- 



solving this first order equation, we obtain ac(o") = exp [a. 



1 



that, for 



Xi 

-X2 

1 

different values of xi and X2 it has a different forms which we divide these forms in the 
following cases: 
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■ By 



V 



I. The case Xi = X2 = 0. 





/ 


1 








\ 




In this case, the curve is in the form acip^ = 







1 







. It is clear that the 




I 




cr 


1 


I 





first column of this matrix is c'((7) and so we have c(cr) = X + (cr , | cr^ , \o^^ for some 
constant K S R^, that its image is analogous to the image of twisted cubic [4J. Also 
the image is similar to the Neil or semi-cubical parabola's graph [7]. The projection of 
this space curve in the direction of z— axis is a parabola in plane. This space curve is 
the simplest curve in R^ under special affine transformations. Its figure is a translation 
of Figured]- (a) by constant K. 
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Theorem 4.1 Space curves with zero special affine curvatures are in the form of 
twisted cubic probably with some translations. 



Figure 1: (a) xi = X2 = 0. (b) xi = 0,X2> 0. 







II. The case Xi = and X2 > 0. 

In this case, we have 

( ' ' 

" ^ sin(^/X2cr) cos(^cj) 

V (cos(0^(7) - 1) ^ sin(^cj) 

So we obtain c{a) = K + (a , cos(^/x2 o") , "^^^^ sin(^/x2 o") + 3^) for K G R^. 
The image of this curve is a translation of Figure [T]-(b) by constant K. Its projection 
in the direction of z— axis is similar to the graph of function cos((t). 



X2 sm(^X2crj 
COs(^/X2Cr) J 



III. The case Xi = and X2 < 0. 

If we use 1x2 1 = — X2 to be the absolute value of X2; in the same way as the previous 
cases, we find that 



/ 



ac{cr)-- 



1 



sm. 



HVbal 



a) 





cosh(/|x2l 



13^{cosh(v/|X2|o-) - 1} 



'\X2\ 



sm 



\ 

(y) v1x2[sinh(v/fx2|o-) 
a) cosh(v/|x2[o-) j 
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Thus c{a) = K+(a, cosh{^/\x2\ cr) , sm'h{^/\x2\ cr) - j^) where K is an 

element of R^. Its image is drown in Figure [2]-(a) probably after a translation. Its 
z— axis projection is similar to the graph of the function cosh((T). 



Figure 2: (a) xi = 0, X2 < 0. (b) xi > 0, X2 = 0. 



IV. The case Xi > and X2 = 0. 

Under these conditions, the adcr) is 



|M+ii? -^xl^^iV^N+M-R) ^xl^^{V^N-M+R) ^ 

j\j^{V3N-M+R) iM + ^R -IxY^iV^N+M-R) 

.-AjsiV^N+M-R) 1 (V3iV-M+i?) iM + ^R 

where 

M = exp(-lxl^'^) cos(^x}^'<^), R = exp(x}/'a) 
iV = exp(-ix}/M smi^xl^'a). 

Therefore with above conditions, we can write 

c{a) = K+{^^{V3N-M+R), {-V^N-M+R), -^(2M+i?) 

for some K G R'^. Its figure is similar to Figure [2]-(b). 
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V. The case Xi < and X2 = 0. 



Such as case III, with use of |x2| = —X2, the conditions of this case lead to the form of 
ac(cj): 

^ iM+^R i|xi|^/^(V3iV+M-i?) l\xl\^/^{^/3N+M-R) \ 

^^(^N+M-R) iM+^R _i|^,|i/3(^Ar_M+i2) 

V ^^{V3N-M+R) ^^(^N+M-R) ^M+^R 
where 

M = exp(i|xir/V) cos(^|xi|^/^/ct), R = exp(-|xir/M 



N 



exp(i|xi|i/V) sin(^|xir/3/ 



And so we have the following curve 

c(a) = /^+(-^(^/3iV+M-i?), -^iyj(^/3iV-M+i?), — ^)2M+ii) 
^3|xi|^/'^ 3|xir/^ 3|xi| 

for some K £ R'^. Its shape is similar to Figure [3]- (a). 



Figure 3: xi <0,X2 = 0. 



VI. The case X11X2 7^ 0. 



In this case, relations are not as simple as previous cases, adcr) in this general case, is 

Bii B12 Bi3 

B21 B22 B23 



in the following form adcr) 



where for 1 < i, j < 3 the entries 



y -B31 -B32 B-33 J 



Bij for brevity are given in Appendix at the end of the paper. 
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Thus the relation c'(cj) = {Bll, B21, B31) signifies the curve c(cj) by integrating of 
the coefficients with respect to a. Therefore we obtain c(cj) = K + (T1,T2,T3) when 
K is an element of R'^ and Tj s are in the forms of variables which are indicated in 
Appendix. 

For different values of constants Xi, X2 7^ 0, there exist various curves and c{a) is 
a translation, contraction or extraction of a curve in the form of these cases. Thus we 
have different figures that each of which is similar to one of the shapes given in Figure 
m (a)-(d). 



Figure 4: (a) Xi,X2> 0. (b) Xi,X2 < 0. (c) Xi < < X2- (d) Xi > > X2- 

Corollary 4.2 In general, every solution o/ Oc : R ^ SL(3, R) is provided by multi- 
plying a special linear transformation and a translation from an acquired curve in above 
cases. In fact, the geometrical sense of above curves can be explained as follows: 
Each curve has two branches. The values of the first and second special affine cur- 
vatures determine "rotation quantities" of the branches that by ascending (descending 
resp.) the values, each branch's bend will increase (decrease resp.). Accordingly, the 
definitions of xi and X2 have geometric interpretations as the usual terminology of cur- 
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vatures. 

In the case of constant xi X2, by using Theorem 3.3, we can classify space 
curves in these cases via special afhne transformations and as a result we have the 
following theorem: 

Theorem 4.3 Each curve of class in satisfied in condition ([7]) with constant 
affine curvatures xi o-nd X2, up to special affine transformations, is the trajectory of a 
one-parameter subgroup of special (unimodular) affine transformations, that is, a curve 
of cases I- VI. 

Finally we give two applications of the classification of space curves by the action 
of special affine transformations and Theorem 4.3: 

Corollary 4.4 In the physical sense, we may assume that each space curve X : 
[a, b] — > is the trajectory of a particle with a specified mass m " in R3 and m the view 
of an observer, that is influenced under the effect of a force F. By the action of special 
affine transformations, particle's path has two conservation laws: (X" x X'") • X"" and 
(X'xX"')-X"", that are, the first and the second special affine curvatures. Therefore, by 
multiplying constant nn? to theses invariants, we find conservation laws as (F x F') • F" 
and (P X F') • F" where P = m • v is the momentum of the particle. If these invariants 
of the trajectory are constant, then the shape of the motion is similar to one of the six 
cases mentioned in the above theorem. 

The derived invariants in above, may have important applications in astronomy, 
fluid mechanics, quantum, general relativity and etc. A reason for this importance is 
that in these areas we deal with the motion of a space particle and it may be of our 
interest to investigate for symmetry properties and invariants of the particle under rigid 
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motions. 



Corollary 4.5 In computer vision and image processing, we may suppose that each 
space curve is one of the characteristic curves on a 3-dimensional object, that are fea- 
sible minimum segment curves that completely signify the object in the viewpoint of 
an observer. Also, if by an effect provided by (orientation-preserving) rotations and 
translations in we change the position of a picture without any change in charac- 
teristic lines, then these cur es will be equivalent under special affine transformation. 
If a characteristic line has constant affine curvatures xi o.^'d X2, then it will be similar 
to one of the cases of curves mentioned in Theorem 4-3. 

For instance, these image invariants provide the most prominent apphcation fields 
in 3D medical imagery, including MRI, ultrasound and CT data, in object recognition, 
symmetry and differential invariant signatures of 3D shapes [HI [9l [12] . 

Appendix 

In case VI of section 4, entries Bij (1 < i,j < 3) are 



^11 



- (72 ci/' X1X2-8 C2 X2 - 144 C2X1 - 192 - 1296 x? 
+8ci^^ X2) + S(l3.85640646ci''^ xl + 12.470765820 xi X2 
+13.85640646 cj-^ C2 X2) + ^( - 8^/^ X2 - 72 ca Xi - 648 x? 
-96 xi - 72 bY^ X1X2-8 cY' C2 X2) } , 



B21 





)} 



- 1.732050808 c 



4/3 



1 
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B22 = 533 = ^{^(l296x? + 144c2Xi+4q/'xi + 36c}/'xiX2 
+4 C2 X2 + 192 X2) + B (62.353831080 Cg^^ xi X2 
+69.28203232 c/^ C2 X2 - 69.28203232 cf^ xl) + ^ (72 C2 Xi 
+96 xi + 648 X? - 36 c'/' xi X2 - M4 cj/' C2 X2 - 4 c'/' xl) } , 

F23 = ^{a(i08cY\1 + 24cl^\l - Gcf\iX2 - 2cY'c2X2 + 12c}/'c2Xi) 
+b(i0.39230485ci/^XiX2 + 3.464101616c}/^C2X2 + 187.0614873 
xcy\l + 41.569219401/^^ + 20.7846097ci/^C2Xi) + d(6cI^\iX2 
-24cJ/'x2 + 2q/'c2X2 - 12q/'c2Xi - 108cJ/'xi)}, 

Which in the above relations we assumed that 

ci = 108x1 + 12(^12xi + 81x?), 
C2 = ^12xi + 81x?, 

A = exp ( - 0.08333333333^^^ — "3^'^^ ^) ^os (0.1443375673^^^ — T/T'^ '^)' 
B = exp{- 0.08333333333^^i — "3^'^^ ^) (0.1443375673^^^ — T/T^ 

2/3 _ , 2 

i:' = exp (0.1666666667^^i 173^ 

A = (cf X2 - 9c}/'^ XI - C2 - 12x^)(cf + 12x2). 

Also Tj s are in the following forms 

Ti = ^{^(5832x?X2 + 864x1X2 - 24cJ/\^_ 81^2/3^3 + 48^2^4 

-24c^/'xixi + 648c2X?X2 - 300^/^02x1 - ^4^^02x1 - 2c^/'c2xi 
-270c}/^X?X2) 41.56921939c}/\i + 140.2961154ci/^x? 

+41.56921939ci/\iX2 - 467.6537182ci''"^x?X2 + 15.58845727ci/^C2X? 
-51.96152424c}/'c2xi + 3.46101616c?/'c2xi) + ^(432xiX2 + 24c2xt 
+81c^/'x? + 2916x?X2 + 24c}/'x2 + 24c?/'xixi + 270c;/'x?X2 
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+324c2xlx2 + 2ci/'c2xi + 9 C2 xl + ^Ocl^"" C2Xixf) } , 
T2 = ^{Ai^864cl^\ix2 + 216c^/' C2X1X2 + l(i2c''/\l + l8cf'c2Xi) 
-5(3367.106771401/^X1X2 - 374.1229746ci/^C2XiX2 
+280.59223086ci/\? + 31.176914544ci/^C2Xi) - i?(l62ci/^x? 
-18c^/^C2Xi - 216^/^02X1X2)}, 
T3 = ^{^(l2ci/\|-324c2X?- 432x1x1 -24c2xi-2916x? + 9c?/' 

•X1X2 + c?/'c2xi + 162cl/'x?X2 + 180^/^02X1X2) + 5(20.784609700^/' 
•X2 - 15.58845727ci/'xiX2 - 1.732050808ci/'c2xi + 280.5922309^/' 
•X?X2 + 31.17691454c}/'c2XiX2) - 5>(l2c}/'x2 - 9c?^^XiX2 - Ci/'c2X^ 
-162c}/'x?X2 - 162C2X? - 216x1x1 - 12c2xi - 1458x? - 18c;/'c2XiX2) }• 

In the latter relations, A, B, D, ci and C2 are the same with relations mentioned in 
above and 

A' = [9c}/\i + q/'c2-cf X2 + 12x^][9c}/'xi + q/'c2 + cf X1 + I2X2] 
x[9c}/'xi + cl/'c2-12x^]. 
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